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ABSTRACT  Direct,  tunable  coupling  between  individually  assembled  graphene  layers  is  a 
next  step  toward  designer  two-dimensional  (2D)  crystal  systems,  with  relevance  for  fundamental 
studies  and  technological  applications.  Here  we  describe  the  fabrication  and  characterization  of 
large-area  (>cm2),  coupled  bilayer  graphene  on  Si02/Si  substrates.  Stacking  two  graphene  films 
leads  to  direct  electronic  interactions  between  layers,  where  the  resulting  film  properties  are 
determined  by  the  local  twist  angle.  Polycrystalline  bilayer  films  have  a  "stained-glass  window" 
appearance  explained  by  the  emergence  of  a  narrow  absorption  band  in  the  visible  spectrum  that 
depends  on  twist  angle.  Direct  measurement  of  layer  orientation  via  electron  diffraction, 
together  with  Raman  and  optical  spectroscopy,  confirms  the  persistence  of  clean  interfaces  over 
large  areas.  Finally,  we  demonstrate  that  interlayer  coupling  can  be  reversibly  turned  off  through 
chemical  modification,  enabling  optical-based  chemical  detection  schemes.  Together,  these  results  suggest  that  2D  crystals  can  be  individually  assembled 
to  form  electronically  coupled  systems  suitable  for  large-scale  applications. 
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The  discovery  and  development  of 
graphene1,2  has  drawn  attention  to 
a  class  of  materials  broadly  defined  as 
two-dimensional  (2D)  crystals.3  As  the  field 
matures  in  understanding  and  exploiting 
single  2D  crystals,  new  research  is  emerging 
to  build  up  materials  from  individual  layers 
forming  "designer"  thin  films.4-7  While  sim¬ 
ple  material  combinations  may  be  realized 
through  direct  growth,8,9  many  others  will 
require  the  physical  stacking  of  individual 
layers.  However,  the  top-down  synthesis  of 
multilayer  structures  is  often  plagued  by  sur¬ 
face  contamination,  which  limits  direct  van  der 
Waals  contact  between  layers.10,11  In  the  sim¬ 
plest  case — the  stacking  of  two  graphene 
layers — new  properties  arise  when  clean  inter¬ 
faces  are  realized.12-14  It  is  already  clear  that 
bilayer  graphene  is  a  diverse  material  system 
because  the  film  properties  vary  as  a  function 
of  relative  orientation  of  one  layer  to  the 
other.  Such  twisted  graphene  systems,15  whose 
theoretical16-18  and  experimental11-14,19-22 
properties  are  now  being  intensely  considered, 
will  guide  our  development  of  other  synthetic 
van  der  Waals  films. 

In  this  work,  we  describe  the  top-down 
synthesis  and  characterization  of  coupled 


bilayer  graphene  films  that  exhibit  intimate 
contact  over  macroscopic  areas  (>cm2). 
The  observation  of  widespread  "colored" 
domains  in  these  stacked  films  reveals  that 
coupling  between  layers  results  in  new 
properties  not  intrinsic  to  the  individual 
components,  an  effect  not  observed  pre¬ 
viously  in  large-area  artificial  bilayer  films  on 
Si02/Si  substrates.  To  determine  if  these 
colored  domains  are  an  intrinsic  feature  of 
the  film,  we  use  Raman  spectroscopy,  which 
has  proven  to  be  useful  in  quantifying 
various  degrees  of  interaction  between  gra¬ 
phene  layers.11,13,19,23,24  Interlayer  hybridiza¬ 
tion  in  twisted  bilayer  graphene  (TBG) 
results  in  measurable  changes  in  the  inten¬ 
sity,  position,  and  shape  of  the  characteristic 
G  (-1600  cm-1)  and  2D  (-2700  cm-1) 
Raman  peaks;  as  such,  a  strong  one-to-one 
correlation  between  Raman  intensity  and 
twist  angle  has  been  established.11,13  In 
addition  to  changes  in  the  Raman  response, 
the  hybridization  of  the  Dirac  cones  in  TBG21 
results  in  changes  of  the  optical  conductiv¬ 
ity,  including  the  emergence  of  an  absorp¬ 
tion  peak11,20  in  the  relatively  wavelength- 
independent  spectrum  of  single-layer  or 
Bernal  bilayer  graphene.25,26  In  this  work, 
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Figure  1 .  Fabrication  and  characterization  of  two  stacked  graphene  films.  (A)  Schematic  of  the  process  flow  used  in  forming 
bilayer  graphene  from  single-layer  CVD  graphene  films.  (B)  Optical  microscope  image  of  a  TBG  film.  The  "red",  "yellow",  and 
"blue"  domains  are  labeled.  (C)  Sample  in  which  the  top  (second)  graphene  layer  is  only  a  partial  layer,  forming  single-layer 
(1 L)  and  bilayer  (2L)  regions.  (D)  AFM  phase  image  (1 5  x  15  /<m2)  showing  the  boundary  between  the  single-layer  and  bilayer 
region  highlighted  by  the  red  box  in  (C). 


optical  spectroscopy  reveals  the  emergence  of  an 
adsorption  peak  in  the  visible  spectrum  for  TBG  do¬ 
mains  with  twist  angles  between  ~10  and  16°,  which 
we  have  independently  confirmed  using  low-energy 
electron  diffraction.  Analogous  to  the  rapid  develop¬ 
ment  of  graphene  research  provided  by  simply  "visua¬ 
lizing"  graphene,27,28  the  visualization  of  specific  twist 
orientations  with  an  optical  microscope  should  further 
enable  the  rapid  study  of  TBG  systems.  As  such,  we  use 
optical  microscopy  to  demonstrate  that  interlayer 
coupling  can  be  effectively  switched  "on"  and  "off" 
through  chemical  functionalization  of  the  top  surface. 

Graphene  films  were  grown  via  low-pressure  chem¬ 
ical  vapor  deposition  (CVD)  in  Cu  foil  enclosures29 
and  transferred  onto  SiO2(100  nm)/Si  substrates 
using  conventional  wet-chemistry  techniques, 30-32 
including  the  use  of  a  poly(methyl  methacrylate) 
(PMMA)  protective  layer33  and  APS  copper  etchant 
100  (Transene  Company,  Inc.).  Before  etching  Cu,  we 
oxygenate  the  Transene  etchant,  which  we  find  results 
in  fewer  carbonaceous  residues  after  PMMA  removal. 
Immediately  after  transferring  the  PMMA/graphene 
film  from  a  H20  bath,  the  substrate  is  spun  at 
2000  rpm  then  4000  rpm  to  remove  bulk  water 
(Figure  1A).  The  sample  is  subsequently  heated  on  a 
hot  plate  (T  =  85  °C,  3  min;  T=  1 50  °C,  20  min)  and  then 
submerged  in  an  acetone  bath,  rinsed  in  acetone  and 
isopropyl  alcohol,  and  dried  with  N2.  Bilayer  samples 
are  generated  by  sequentially  transferring  a  second 
CVD  graphene  layer  onto  the  graphene/Si02/Si  sub¬ 
strate  and  repeating  the  steps  described  above.  The 
polycrystalline  nature  of  the  starting  CVD  graphene 
results  in  bilayers  with  regions  of  varying  twist  angle 
and  is  referred  to  here  as  twisted  bilayer  graphene. 

RESULTS  AND  DISCUSSION 

Optical  inspection  of  as-fabricated  TBG  samples  re¬ 
veals  expected  features  such  as  wrinkles  and  folds,34,35 
together  with  an  unexpected  patchwork  of  colored 


domains  that  vary  in  size  and  shape  across  the  sample. 
Figure  1 B  shows  an  optical  microscope  image  (Olympus 
DP25  CCD  camera)  of  a  TBG  film  with  distinct  regions 
that  appear  red,  yellow,  and  blue.  No  filters  are  neces¬ 
sary  to  see  these  colored  features  (see  Supporting 
Information).  Notably,  subsequent  thermal  annealing 
(e.g.,  400  °C  in  Ar/H2)36  or  additional  solvent  cleaning 
does  not  change  the  extent  of  these  colored  regions, 
suggesting  that  their  presence  is  not  related  to  extrin¬ 
sic  processing  residues.  Figure  1C  shows  a  higher 
resolution  image  of  a  different  TBG  sample  in  which 
the  second  (top)  graphene  layer  is  noncontinuous. 
Flere  we  find  that  the  colored  domains  are  only 
observed  in  the  bilayer  regions  and  not  the  neighbor¬ 
ing  single-layer  regions.  In  addition  to  these  macro¬ 
scopic  colored  features,  small  isolated  islands  are 
observed  (bluish  "dots"  Figure  1C;  discussed  below), 
together  with  the  wrinkles/folds  mentioned  earlier. 

Atomic  force  microscope  (AFM)  imaging  at  the 
boundary  between  a  bilayer  and  single-layer  region 
shows  an  increased  level  of  roughness  in  the  bilayer 
region  (Figure  IE).  This  increased  roughness  is  in  the 
form  of  small  isolated  islands  that  have  typical  foot¬ 
prints  oflessthan  1  /im2and  heightof ~10— 40 nm and 
are  separated  by  relativity  smooth  areas.  Similar  islands 
have  been  observed  elsewhere  when  multiple  gra¬ 
phene  layers  are  individually  assembled,10'13  where 
the  largest  islands  align  well  with  the  bluish  dots 
observed  optically  (Figure  1C).  Cross-sectional  imag¬ 
ing  of  stacked  graphene  layers  confirms  the  pres¬ 
ence  of  locally  trapped  amorphous  material  in  such 
samples,  separated  by  pristine  regions  in  direct  van  der 
Waals  contact.10  The  ability  of  trapped  (interlayer) 
hydrocarbons/adsorbates  to  segregate  over  micro- 
meter-scale  distances  gives  rise  to  these  isolated 
islands.10  We  note  that  using  a  similar  sample  fabrica¬ 
tion  procedure  as  described  above,  we  have  measured 
coupling  between  atomically  flat  epitaxial  graphene 
on  SiC  and  transferred  CVD  graphene  films,21'37 
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supporting  the  presence  of  clean  interface  regions 
between  such  isolated  islands. 

The  direct  correlation  between  the  Raman  spectral 
response  and  layer  orientation  in  TBG  films11'13  pro¬ 
vides  a  rapid  means  to  determine  if  interlayer  coupling 
is  present,  as  well  as  to  coarsely  identify  twist  angle  in 
coupled  bilayers.  Arguably  the  most  dramatic  feature 
in  TBG  Raman  spectra  is  a  substantial  enhancement 
(>20x)  in  the  G  peak  intensity  when  the  laser  photon 
energy  (£photon)  Is  approximately  equal  to  the  energy 
level  where  the  Dirac  cones  overlap/hybridize.  Because 
the  position  of  band  overlap  is  dependent  on  twist 
angle  (9),  the  orientation  at  which  this  maximum 
Raman  enhancement  occurs  is  dubbed  the  critical 
angle  (0C),  with  a  critical  energy  E0c.  Conveniently,  this 
G  peak  enhancement  occurs  in  a  somewhat  narrow 
range — when  6  is  within  a  few  degrees  (about  ±2°)  of 
0C  or  when  £photon  is  within  a  few  hundred  meV  of  Eoc- 
When  outside  of  this  window,  the  Raman  spectra  can 
only  indicate  if  9  <  0C  or  9  >  9C.UA3 

Unlike  two  noninteracting  stacked  graphene  sheets, 
the  Raman  spectra  for  our  TBG  samples  show  a  rich 
variation  in  peak  intensities  and  shapes  as  would  be 
expected  for  two  interacting  layers  with  distinct  twist 
angles.  Figure  2  shows  Raman  maps  and  spectra  from 
TBG  films  in  Figure  1,  which  contain  red,  yellow,  and 
blue  domains.  Figure  2A,B  shows  a  map  of  the  G/2D 
peak  ratio  (integrated  intensity)  for  the  TBG  film  in 
Figure  1 C,  measured  at  two  different  wavelengths  (488 
and  532  nm).  Close  inspection  of  the  map  reveals  six 
unique  intensity  ratios  in  the  TBG  region  (labeled). 
Notably,  the  colored  domains  in  the  optical  image 
align  very  well  with  the  domains  imaged  by  Raman. 
This  result,  together  with  many  other  Raman/optical 
image  comparisons  (not  shown),  provides  strong  sup¬ 
port  for  the  colored  domains  being  related  to  an 
intrinsic  property  of  the  film  and  being  uniquely 
dependent  on  the  twist  angle  between  the  graphene 
layers. 

Individual  Raman  spectra  (Figure  2C,D)  from  the 
samples  shown  in  Figure  1B,C  can  be  understood  in 
the  frameworks  recently  presented  by  Flavener  et  al.u 
and  Kim  et  al ,13  The  use  of  two  different  photon 
energies  allows  us  to  order  the  spectra  in  terms  of 
increasing  twist  angle  through  analyzing  the  G/2D 
peak  ratios  and  peak  positions.  A  strong  G  peak 
enhancement  occurs  at  the  yellow  domain  using 
fphoton  =  2.54  eV  (488  nm)  and  at  the  red  domain  using 
fphoton  =  2.33  eV  (532  nm).  If  we  assume  £photon  here 
is  equal  to  E0c,  this  corresponds  to  0yeiiow  14.5°  and 
0red  ps  12.50.11  Comparing  histograms  of  the  normal¬ 
ized  integrated  G  peak  intensity  (AG)  for  many  spectra 
within  these  same  regions  sheds  light  on  the  closeness 
of  £photon  to  E0c,  as  well  as  on  the  angular  distribution  or 
deviation  in  9  (Figure  2E,F).  For  example,  at  the  reso¬ 
nance  condition  (£photon  =  E0c)  for  9C=  1 2.5°,  the  inten¬ 
sity  of  Ag  decreases  by  50%  within  approximately  ±1° 
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Figure  2.  Raman  mapping  and  spectra  of  TBG  domains 
using  a  (A,C)  488  nm  (2.54  eV)  and  (B,D)  532  nm  (2.33  eV) 
laser.  The  map  is  of  the  TBG  region  shown  in  Figure  1 C.  (C,D) 
Averaged  spectra  from  different  regions  of  the  sample  as 
labeled.  The  same  areas  are  analyzed  for  both  wavelengths. 
The  spectra  are  offset  for  clarity  and  are  arranged  in  order  of 
increasing  twist  angle.  The  G  and  2D  peaks  are  labeled.  (E,F) 
Flistogram  of  the  normalized  integrated  G  peak  intensity 
(AG)  taken  from  individual  spectra  within  the  yellow  and  red 
regions  shown  in  (A)  and  (B)  at  each  wavelength  (labeled). 
(E)  For  the  red  distribution:  mean  =  46.7,  median  =  48.1, 
range  =  46.2.  (F)  For  the  yellow  distribution:  mean  =  18.6, 
median  =  17.8,  range  =  37.3. 


and  by  80%  within  approximately  ±2.5°  twist  varia¬ 
tion.11  If  the  twist  angle  of  the  red  domain  equals  12.5° 
with  a  small  twist  deviation  (±1°),  the  histogram  at  2  = 
532  nm  should  appear  like  that  in  Figure  2E,  having  a 
high  mean  value  and  left  skew.  On  the  other  hand,  if  9 
(fphoton)  differs  by  even  1 0  (1 50  meV)  from  9C  (E0c),  then 
the  Aq  intensity  distribution  will  have  a  lower  mean 
value  (by  up  to  ~50%)  and  have  a  more  symmetric  or 
right  skew  as  observed  with  the  yellow  domain  for  2  = 
488  nm  (Figure 2F). Together,  this  implies  9redm  12.5  ±  1° 
and  that  0yeiiOw  is  slightly  larger  than  1 4.5°.  Finally,  for  the 
blue  domain,  the  slight  increase  in  the  G/2D  ratio  and 
upshift  of  the  2D  peak  from  2  =  488  to  532  nm  indicates  it 
is  by  several  degrees  smaller  than  0red. 

The  presence  of  colored  TBG  domains  suggests 
there  are  unique  differences  in  the  absorption  spec¬ 
trum  that  are  dependent  on  twist  angle.  The  electronic 
hybridization  in  TBG  leads  to  changes  in  the  optical 
transition  matrix  elements  and  joint  density  of  states 
(JDOS),  which  can  strongly  affect  the  absorption  prop¬ 
erties.1  1,20  Figure  3  shows  the  measured  and  calculated 
optical  contrast  spectra  for  a  red,  yellow,  blue,  and 
large-angle  bilayer  domain  (9  >  9C  for  £photon  =  2.54  eV), 
together  with  a  single-layer  region.  The  contrast  spec¬ 
trum  is  defined  as  C(2)  =  [fi0(2)  —  R(2)]//?0(2),  where 
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Figure  3.  Averaged  contrast  spectra  (C(2))forTBG  domains. 
(A)  Measured  C(2)  for  a  red,  yellow,  blue,  and  large-angle 
(0  >  0C)  domain,  together  with  a  single-layer  region.  All 
spectra  are  referenced  to  the  Si02(1 00  nm)/Si  substrate.  The 
calculated  contrast  for  single-layer  (1L)  and  two-layer  (2L) 
graphene  is  included.  The  red,  yellow,  and  0  >  0C  domain 
were  measured  from  Figure  1C.  (B)  Contrast  difference  of 
the  yellow,  red,  and  blue  domains  with  respect  to  the  large- 
angle  domain,  together  with  the  best-fit  curve  calculated 
using  £2(<o)  shown  in  the  inset. 

R0[X)  is  the  reflection  spectra  from  the  substrate  (Si02/ 
Si)  and  R(X)  is  reflection  from  the  film  plus  substrate28,38 
(see  Supporting  Information).  For  single-layer  gra¬ 
phene,  we  find  a  reasonably  good  match  between 
the  measured  and  calculated  C(2)  using  an  index  of 
refraction  for  graphene  of  ng  =  2.6  —  1.3/28  and  thick¬ 
ness  of  0.34  nm,  together  with  literature  values  of  n(X) 
for  Si02  and  Si.39  For  two-layer  graphene,  the  mea¬ 
sured  C(l)  (TBG  domains)  and  calculated  C(A)  (Bernal 
stacking)  do  not  agree  well.  In  particular,  C(2)  for  the 
colored  TBG  domains  is  not  the  smooth  envelope 
shape,  but  instead  there  is  a  shoulder  feature  unique 
to  each  domain. 

Since  contrast  here  is  approximately  proportional  to 
absorption,11,20  the  difference  between  two  contrast 
spectra  highlights  relative  absorption  features  of  one 
domain  to  another.  Flence  we  define  a  contrast  difference 
spectrum,  dQ(/l)  =  C,(l)  —  Ce>gc(X),  where  the  subscript 
i  indicates  yellow,  red,  or  blue.  We  chose  a  large-angle 


domain  as  a  common  reference  since  its  properties  in 
visible  wavelengths  are  close  to  that  of  two  independent 
graphene  layers.  Figure  3B  compares  dC(X)  among  the 
red,  yellow,  and  blue  domains.  For  each  curve,  there  is  a 
distinct  peak  (Epeak),  indicating  enhanced  relative  absorp¬ 
tion  at  specific  regions  of  the  visible  spectrum. 

Given  knowledge  of  the  substrate  (see  Supporting 
Information),  it  should,  in  principle,  be  possible  to 
model  the  complex  dielectric  function  e(a>)  =  en (a>)  — 
i£2(ru)  of  TBG  that  produces  the  strong  features  in  the 
contrast  difference  data  (e2(a>)  is  directly  related  to  the 
intrinsic  material  absorption  and  includes  contribu¬ 
tions  from  the  TBG  hybridization  states).  To  carry  out 
the  modeling,  we  used  a  commercial  software  package 
able  to  model  the  reflectance  from  layered  stacks.40  We 
were  able  to  obtain  rough  agreement  by  using  the 
time-honored  approach  of  modeling  narrow  absorp¬ 
tions  with  Gaussian  or  Lorentzian  oscillators  but  found 
that  the  fits  improved  significantly  when  we  employed 
a  more  flexible  model  that,  in  a  Kramers— Kronig  con¬ 
sistent  fashion,  reproduces  the  van  Fiove  singularity 
(vFHs)  and  linear  density  of  states  expected  for  the  TBG 
absorption.11,41  Our  fits  to  dC{X)  for  each  domain  angle 
are  shown  with  dashed/dotted  lines  in  Figure  3B,  and 
the  corresponding  e2(X)  values  are  included  in  the  inset 
of  Figure  3B.  While  qualitative  in  nature,  the  fit  repro¬ 
duces  the  small  dip  after  the  main  peak,  as  well  as  the 
increasing  tail  at  lower  wavelengths,  indicating  that 
the  singularity  plus  linear  shape  of  the  JDOS  is  a  good 
approximation  for  the  presence  and  shape  of  the 
absorption  feature.  In  addition,  the  fwhm  of  the  ab¬ 
sorption  feature  is  approximately  0.25  eV,  close  to  that 
predicated  by  DFT  calculations.11  Assuming  here  that 
£peak  equals  ESc,  this  corresponds  to  ©b|ue  «  10.5°,  0rec i 
~  1 2.5°,  0yeiiow  ^  15.50,11  in  excellent  agreement  with 
the  Raman  results  and  analysis  in  Figure  2. 

To  directly  quantify  interlayer  twist  angles  and  var¬ 
iations  within  one  layer's  in-plane  orientation,  we 
conduct  microdiffraction  experiments  using  low-en¬ 
ergy  electron  microscopy/diffraction  (LEEM/LEED).  Un¬ 
like  LEED  of  TBG  on  SiC,37  the  LEED  spots  of  graphene 
and  TBG  on  Si02  are  diffuse  (Figure  4A,B)  due  to  the 
nanometer-scale  substrate  roughness  of  Si02.42,43 
Since  diffracting  electrons  from  the  bottom  layer  in 
TBG  on  Si02  are  greatly  attenuated  in  these  samples, 
we  presume  that  the  orientation  of  the  diffraction 
pattern  represents  primarily  that  of  the  top  surface 
layer.  As  such,  we  only  measure  one  family  of  diffrac¬ 
tion  spots  on  either  TBG  or  single-layer  graphene 
(Figure  4A,B). 

Mapping  diffraction  patterns  allows  direct  correla¬ 
tion  between  the  angular  orientation  and  spatial  ex¬ 
tent  of  TBG  domains.  Figure  4C  shows  such  a  false  color 
map  of  LEED  angular  orientation.  The  measurement 
was  carried  out  by  translating  the  partial  TBG  film 
(shown  in  Figure  4D)  with  respect  to  the  electron  beam 
(5  //m  diameter)  over  250  x  200  firr\2  (50  x  40  data 
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Figure  4.  LEED  analysis  of  TBG  films.  (A)  Typical  LEED  pattern  of  a  TBG  film  measured  using  5  /im  diameter  area  selective 
aperture  and  an  electron  energy  of  55  eV.  (B)  Intensity  profile  along  the  dotted  arc  in  (A).  (C)  False  color  map  of  the  LEED 
angular  orientations  of  a  TBG  film.  (D)  Optical  microscope  images  of  the  TBG  film  measured  in  (C).  The  numbers  indicate  the 
approximate  LEED  orientations.  (E)  Vertical  line  scans  taken  at  the  black  dash  line  in  (C).  (F)  Histogram  from  the  angular 
contour  plot  in  (C).  The  red  and  neighboring  single-layer  domain  from  (C)  are  specified  in  (F). 


points)  .  We  specify  the  orientation  of  the  diffraction 
pattern  from  the  angle  defined  in  Figure  4A,  based  on 
the  intensity  profile  along  the  dotted  arc  (Figure  4B). 
Distinct  angle-specific  domains  (labeled)  in  both  the 
single-layer  and  bilayer  regions  become  apparent  as 
highlighted  by  the  black  outlines  in  Figure  4C,D.  Cross- 
sectional  line  scans  across  TBG  and  the  neighbor¬ 
ing  single-layer  regions  allows  us  to  quantify  the 
relative  twist  angle  between  the  upper  and  lower 
layers  (Figure  4E).  A  histogram  of  all  diffraction  spots 
also  reveals  the  relative  LEED  orientation,  as  well  as  the 
deviation  in  graphene's  in-plane  orientation  within  a 
domain.  Each  peak  (or  cluster)  in  Figure  4F  represents  a 
domain  within  Figure  4C  and  shows  a  typical  domain 
angular  variation  of  ±1°.  Using  this  same  process  to 
determine  twist  angles  of  other  colored  domains  (see 
Supporting  Information),  we  find  0b|ue  =11  ±  1  °,  0red  = 
1 3  ±  1  °,  and  0yeiiow  =  1 5  ±  1  °,  in  good  agreement  with 
Raman  and  optical  spectroscopy  measurements  de¬ 
spite  the  large  angular  uncertainties  due  to  broadened 
diffraction  patterns. 

The  macroscopic  appearance  of  primarily  three  dis¬ 
tinct,  colored  TBG  domains  is  striking  since  the  energy 
of  the  vhls  feature  continuously  varies  with  twist 
angle.12  As  such,  this  may  be  an  indication  that  the 
fabrication  process  leads  to  preferred  TBG  orientations. 
A  larger  optical  microscope  survey  does  show  that 
some  TBG  domains  appear  to  have  distinct,  mixed 
color  phases  (e.g.,  Figure  IB  and  Figure  5B).  We  are 
currently  conducting  higher  spatial  and  spectral  reso¬ 
lution  imaging  to  determine  if  there  are  statistically 
significant  differences  in  the  areal  coverage  of  specific 
TBG  domains.  It  can  be  said  that  the  observed  mixed 
colored  domains  are  at  least  in  part  due  to  the  mea¬ 
sured  angular  (or  twist)  deviation  within  each  domain, 
where  a  few  degree  in-plane  rotations  can  result  in  a 
200—300  meV  shift  in  the  absorption  feature  (Figure  3), 
as  well  as  E0c.  As  such,  it  is  most  common  to  find  mixing 
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Figure  5.  Local  addition  (A)  and  removal  (B)  of  fluorine 
on  the  top  surface  of  TBG  films.  Lines  help  guide  the  eye 
where  TBG  was  fluorinated  (labeled  "FG")  and  reduced 
(labeled  "rFG"). 

of  blue  (Sbiue  =  1 1  ±  1  °)  and  red  (0red  =  1 3  ±  1  °),  or  red 
(0red  =  1 3  ±  1  °)  and  yellow  (0yeiiow  =  1 5  ±  1  °)  domains, 
but  not  the  blue  (0b|ue  =11  ±  1  °)  and  yellow  (0yeiiow  = 
15  ±  1°)  domains,  which  have  an  angular  separation 
larger  than  the  angular  deviation  of  domains  in  these 
samples.  Overall  improvements  in  TBG  domain  size  could 
be  greatly  facilitated  through  the  use  of  highly  orientated 
graphene  films  as  recently  demonstrated  on  Au-foil 
substrates,44  while  local  twist  angle  variations  may  be 
reduced  through  further  optimizing  the  transfer  process, 
including  the  use  of  smoother  starting  surfaces  and  low 
surface  energy  solvents  during  transfer  and  drying.45 

Finally,  we  demonstrate  lithographic  control  of  in¬ 
terlayer  coupling  by  selective  fluorination  of  the  top 
graphene  layer  using  XeF2  gas.46  Following  fluorina¬ 
tion,  even  for  very  short  XeF2  exposures  (e.g.,  1  s,  1  Torr 
XeF2),  color  contrast  is  fully  quenched  in  TBG  domains 
with  twist  angles  between  10  and  16°.  This  is  seen  in 
Figure  5A  following  fluorination  in  a  pattern  defined  by 
a  PMMA  mask.  Upon  fluorine  desorption  (7=  1 75  °C  for 
1  h  in  flowing  Ar),  the  color  contrast  returns,  demon¬ 
strating  the  sensitivity  of  interlayer  coupling  to  cova¬ 
lent  functionalization  and  the  usefulness  of  optical 
characterization  as  a  probe  for  interlayer  hybridization. 
The  reappearance  of  interlayer  coupling  after  defluo¬ 
rination  is  noteworthy  since  defects  are  often  intro¬ 
duced  during  the  adsorption/desorption  process,46 
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highlighting  the  overall  robust  nature  of  coupling. 
Fluorine  adsorption  changes  the  electronic  properties 
of  graphene  by  reducing  the  charge  in  the  conducting 
n  orbitals,  by  introducing  scattering  centers,  and  by 
opening  band  gaps.46  Fluorine  also  structurally  de¬ 
forms  the  graphene  skeleton  as  bond  angles  shift 
due  to  transiting  sp2  to  sp3  carbon  bonding,  which  is 
confirmed  by  the  emergence  of  the  D  peak  in  Raman 
spectroscopy46,47  (not  shown).  Thus,  decoupling  via 
mechanical  separation  or  via  quenching  of  the  hybri¬ 
dized  electronic  state  are  promising  routes  for  optical- 
based  chemical  detection  using  TBG  films  with  twist 
angles  between  10  and  16°.  In  addition  to  chemical  or 
mechanical  routes,  it  is  likely  that  electrostatic  gating 
will  provide  tunable  control  over  interlayer  coupling 
and  band  structure,  similar  to  that  reported  for  Bernal 
(AB-stacked)  bilayer  graphene.48,49 

CONCLUSION 

In  summary,  we  show  that  interlayer  coupling  be¬ 
tween  two  large-area,  individually  stacked  graphene 


films  is  possible  and  that  interlayer  coupling  enables 
the  direct  visualization  of  specific  bilayer  orienta¬ 
tions  on  Si02/Si  substrates.  The  diffusion  of  interlayer 
adsorbates  into  localized  islands  facilitates  the  forma¬ 
tion  of  atomically  clean  regions  that  can  electronically 
hybridize.  The  polycrystalline  nature  of  the  initial  CVD 
graphene  leads  to  bilayer  films  whose  properties  de¬ 
pend  on  the  relative  twist  angle  between  layers.  The 
interlayer  coupling  results  in  the  emergence  of  a 
relatively  narrow  absorption  peak  in  the  visible  range 
for  twist  angles  between  10  and  16°,  where  the  twist 
angle  variations  within  a  single  domain  are  found  to  be 
approximately  ±1°  in  these  samples.  By  selectively 
functionalizing  the  top  surface  of  these  bilayer  films, 
the  interlayer  coupling  is  effectively  quenched,  open¬ 
ing  possibilities  of  optical-based  chemical  detection 
applications.  We  believe  the  results  presented  here  will 
facilitate  the  formation  of  various  large-area  homo- 
and  heterogeneous  bilayer  systems,  where  interlayer 
hybridization  can  result  in  exploitable  electronic  or 
optical  properties. 


METHODS 

Raman  measurements  were  performed  using  a  confocal 
geometry.  Dichroic  beam  splitters  were  used  to  reflect  single¬ 
mode  488  or  532  nm  laser  light  onto  the  excitation/detection 
optical  axis.  A  lOOx  microscope  objective  (NA  =  0.65)  focused 
the  laser  (spot  ^0.4  /<m)  onto  the  sample  and  gathered  Raman 
scattered  light  for  detection.  The  Raman  scattered  photons 
were  dispersed  in  a  half-meter  Acton  Sp-2500  spectrometer 
and  were  detected  using  a  Princeton  Instruments  CCD  array 
(Spec-1 0:400BR  back-thinned,  deep-depleted  array). 

Optical  spectroscopy  measurements  were  performed  with  an 
inverted  microsope  (Nikon  TE2000)  coupled  directly  through  a 
side  port  to  an  imaging  spectrometer  (Princeton  Instruments 
MicroSpec  300).  The  spectrometer  was  equipped  with  an  8 /< m 
pixel  CCD  (Andor  model  885)  and  entrance  slit  assembly  that 
could  be  temporarily  moved  aside  for  imaging.  Spatial  registra¬ 
tion  between  images  and  spectra  was  assured  because  the 
spectrometer  grating  (50  groove/mm)  could  be  exchanged 
under  computer  control  with  a  mirror  for  direct  imaging.  The 
sample  was  illuminated  through  the  20x  objective  (NA  =  0.45) 
with  a  quartz-tungsten-halogen  lamp.  The  Koehler  illuminator 
was  apertured  to  restrict  the  incident  rays  to  be  within  ±10°  of 
the  sample  normal.  Reflected  light  was  collected  by  the  same 
objective.  For  optical  modeling  of  the  contrast  difference 
spectra,  to  reproduce  the  expected  van  Hove  singularity  and 
linear  density  of  states  in  TBG,  we  used  the  software's40  built-in 
functions  psemiO  and  psemil. 40,50,51 

To  slightly  improve  the  diffraction  contrast  during  LEED 
measurements,  we  first  deposited  a  self-assembled  mono- 
layer  of  hexamethyldisilazane  (HMDS)  on  Si0252  before  TBG 
deposition,  in  order  to  reduce  surface  roughness.  We  then 
acquired  the  LEED  images  via  a  5  /urn  illumination  aperture 
selecting  the  single  graphene  layer  and  then  the  TBG  do¬ 
mains  (see  Supporting  Information).  In  order  to  ensure  that 
both  single-layer  and  bilayer  regions  have  uniform  orienta¬ 
tions  in  the  field  of  view,  dark-field  LEEM  images  were 
acquired,  as  shown  in  Supporting  Information  Figure  S4. 
Fluorination  experiments  were  carried  out  using  a  Xactix 
XeF2  etching  tool.  In  Figure  5A,  the  sample  was  exposed  in 
pulse  mode  with  the  following  parameters:  10  cycles,  30  s/ 
cycle,  1  Torr  XeF2,  35  Torr  N2. 
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